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Abstract 

Quantum mechanics has been formulated in phase space, with the 
Wigner function as the representative of the quantum density opera- 
tor, and classical mechanics has been formulated in Hilbert space, with 
the Groenewold operator as the representative of the classical Liouville 
density function. Semiclassical approximations to the quantum evolu- 
tion of the Wigner function have been defined, enabling the quantum 
evolution to be approached from a classical starting point. Now anal- 
ogous semiquantum approximations to the classical evolution of the 
Groenewold operator are defined, enabling the classical evolution to 
be approached from a quantum starting point. Simple nonlinear sys- 
tems with one degree of freedom are considered, whose Hamiltonians 
are polynomials in the Hamiltonian of the simple harmonic oscilla- 
tor. The behaviour of expectation values of simple observables and 
of eigenvalues of the Groenewold operator, are calculated numerically 
and compared for the various semiclassical and semiquantum approx- 
imations. 
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1 Introduction 



Developments in nanotechnology have focussed increasing attention on the 
interface between quantum mechanics and classical mechanics. Since the 
earliest days of the quantum theory, the nature of this interface has been 
explored by means of various semiclassical approximations with the intro- 
duction of quantum corrections to classical behaviour, characterised by ex- 
pansions of quantities of interest in asymptotic series in increasing powers 
of Planck's constant h ; see for example [U El El E] • One such approach 
focusses on Wigner's quasiprobability density, which is the representative 
of the quantum density operator in the formulation of quantum mechanics 
on phase space, and considers quantum corrections to the classical Liouvil- 
lean evolution. This approach has its roots in Wigner's paper £Q where the 
quasiprobability function was first introduced, and it has recently been de- 
veloped in some theoretical detail by Osborn and Molzahn |3J, working in 
the Heisenberg picture. 

Recently an alternative approach has been suggested, which we refer to 
as semiquantum mechanics jH]. This considers the quantum-classical inter- 
face from the other side, that is to say, from the quantum side. The idea 
is to consider a series of classical corrections to quantum mechanics, again 
characterised by increasing powers of h. To facilitate this, we first refor- 
mulate classical mechanics, representing classical observables on phase space 
by hermitian operators on Hilbert space [Hj, and in particular representing 
the Liouville density by a quasidensity operator [3 El El (the Groenewold 
operator). This has many of the properties of a true density operator but 
is not positive-definite in general It is the analogue in the Hilbert space 
formulation of classical mechanics, of the Wigner function in the phase space 
formulation of quantum mechanics. 

The evolution in time of the Groenewold operator corresponds to the 
classical evolution of the Liouville density, and has been given in [Hj as a 
series of terms in increasing powers of H, with the first term corresponding to 
the familiar quantum evolution. It is consideration of these terms successively 
that defines a series of approximations to classical dynamics, starting with 
quantum dynamics at the lowest order, and hence defining what we mean by 
semiquantum mechanics. 

In what follows, we shall apply these ideas to some very simple non- 
linear systems in one dimension. These are necessarily integrable, and so 
incapable of showing interesting dynamical behaviour such as chaos. How- 
ever the nonlinearity provides a preliminary testing ground for semiquantum 
mechanics, as we consider successive classical corrections to the quantum 
dynamics, and see how the classical behaviour of expectation values of key 



observables emerges. We are able to identify some differences between semi- 
quantum and semiclassical approximations at the interface between quantum 
and classical behaviours. 

2 Semiquantum mechanics 

The formulation of quantum mechanics in phase space is determined with 
the help of the unitary Weyl-Wigner transform W, which maps each quan- 
tum observable (hermitian operator) A on Hilbert space into a real-valued 
function A on phase space. In the case of one degree of freedom, and when 
A is regarded as an integral operator with kernel Ak{x,v) in the coordinate 
representation, the action of W is defined by ITU IT2"] 

A(q, p) = W(A){q,p) = J A K (q-x/2, q + x/2)e ipx/h dx. (1) 

In the particular case of the density operator p, we obtain the Wigner function 

W(q,p)=W(p/2irh)(q,p). (2) 

Conversely, the formulation of classical mechanics in Hilbert space has 
been defined jHJ E] using the inverse transform W _1 , which is Weyl's quan- 
tization map. This acts on a classical observable (a real-valued function on 
phase space) A to produce an hermitian operator A = W~ 1 (A) with kernel 
A K given by 

A K (x,y)=W~ 1 (A) K (x,y) = 7 ±- [ A([x + y]/2, p) e^~^ h dp . (3) 

2irn J 

Typically W _1 is used only to establish the quantization of a given classical 
system. However, we can if we wish use it to map all of classical mechanics, 
including classical dynamics, into a Hilbert space formulation [HIE]. Then h 
in (J3J) can be thought of as an arbitrary constant with dimensions of action 
jH], to be equated with Planck's constant if and when desired. 

Given a Liouville probability density on phase space, that is to say a 
function p(q, p) satisfying 

J P{Q, p)dqdp=l, p(q, p) > , (4) 

we define the corresponding hermitian (integral) operator G = W~ 1 (27ihp) 
using (JHJ). This operator G is the Groenewold operator, which has all the 



properties of a true density operator except that it is not in general positive- 
definite. Thus if A = W~ l {A), then 

Tr(G) = l, Tr(G 2 )<l, 



(A) = J p(q, p) A(q, p) dqdp = (A) = Tr(Gi) , (5) 

but not all eigenvalues of G need be positive. 

It was shown in jHj that for suitably smooth Hamiltonians, the classical 
time evolution of the Liouville density, 

dp/dt = p q H p - p p Hq (6) 

is mapped into the evolution 

Q t - = ^ ~ 24 ^P'^ ~ ^H qp , G qp ] + [H pp , G qq 



([Hqqqq, G ppp p\ 4[H qqqp , G qppp ] + 6 [H qqpp , G qqpp ] 

^[H qppp , G qqqp \ -\- [H p ppp, G qqqq \ ) -\- . . . (7) 



where the classical Hamiltonian H(q, p) is represented by the operator H = 
W^ 1 (H) on Hilbert space. In (jHJ) and (J2J) we have introduced the notation 

A p = dA/dp , A q = dA/dq , A qp = d 2 /dqdp , 



A p = w-\a p ) = [q, A]/(ih) , A q = w-\A q ) = [A,p]/(iH) , 

Aqp = A pq = W~ l (A qp ) = i{q,A],p]/(ih)\ (8) 

and so on. In ( |7| ). the numerical coefficients are those in the expansion of 
l/(2ism(h/2)) in ascending powers of h. See jH] for details. 

Thus, to lowest order in h, the classical evolution of G is given by the 
quantum equation 



(9) 



whereas the first order semiquantum approximation has 

Q t - = ^[#> G \ - 24 ([^ffli G pp ] - 2[H qp , G qp ] + [H pp , G qq ]j , (10) 
and so on. 

In what follows, for simple nonlinear Hamiltonians H, we consider first 
and second order semiquantum approximations to the classical evolution of 
G, and look at the resulting effects on expectation values of key observables 
by using those approximations to G in (J^J). We also consider the behaviour 
of the spectrum of G as these successive approximations are introduced. 

Finally, we compare these results with corresponding results obtained 
using semiclassical approximations that are determined by considering suc- 
cesively more terms from the well-known series analogous to (j2J) for the evo- 
lution of the Wigner function W(q, p, t), namely ^3] 

dW 

{H, W}* = (H g W p - H p W q ) 



Of 

h 2 



3!2 2 



(HqqqW ppp 3HqqpWq pp ~\~ <5Hq pp Wqq p H ppp W qq q^ 



5!2 4 



ri«/ ( H qqqqqW ppppp 5Hqqqq p Wq pppp ~\~ iOHqqq pp Wq qppp ^-OHqqppp^Vqqqpp 



~\~^Hq ppp ^Wqqqq p HpppppW^qqqqq) ■ ■ ■ , (H) 

where { , }* denotes the Moyal (star) bracket. The numerical coefficients in 
this series are those in the expansion of 2 sin (ft/2) /ft in ascending powers of 
ft. 

In this case, the formulas for the quantum averages, analogous to the 
classical formulas (0), are 



(A)(t) = Tr (p(t)i) = (A)(t) = / W(q,p,t)A(q,p)dqdp, (12) 



where p is the quantum density operator, and A = W(A). Semiclassical 
approximations to these averages are obtained by inserting in (fT2j) . successive 
approximations to W found from (JTTJ). 



3 A class of simple nonlinear systems 

The equations of motion for linear systems (quadratic Hamitonians) in the 
quantum and classical regimes are identical, whether represented in the 



Hilbert space or phase space formulations. Differences only arise for non- 
linear systems which, in most cases, can only be studied numerically. 

One class of nonlinear dynamical systems for which some analytic results 
can be obtained are systems for which the Hamiltonian is a polynomial in 
the simple harmonic oscillator Hamiltonian 

H = p 2 /2m + muj 2 q 2 /2. (13) 

We introduce the set of classical Hamiltonians of the form 

K 

H = E^bk (H /E) h , (14) 

fc=0 

where the bk are dimensionless constants, K is a positive integer and E is a 
suitable constant with the dimensions of Energy. Hamiltonians of this form 
may seem somewhat artificial, but interactions of this general form play a 
significant role in the theory and simulation of Kerr media and laser-trapped 
Bose-Einstein condensates |14| I15j. 

The primary advantages of these Hamiltonians H and their quantizations 
H = W -1 (if) in the present context are that they generate exactly solvable 
dynamics in both cases, and that the classical dynamics and the quantum 
dynamics display very different behaviours for K > 1, for sufficiently large 
times [TH] . 

In order to investigate the classical dynamics generated on phase space by 
such Hamiltonians, we first introduce the dimensionless complex conjugate 
variables 

a = (y/muq + ip / \/muj) / V2~R , a = (y/muq — ip / 'y 'muj) j 'V2h , (15) 

so that H = hwaa. In the classical context, h should be considered as 
simply an arbitrary constant with dimensions of action. However, we think 
of E as a 'classical' energy, and fkv as a 'quantum' energy, and define the 
dimensionless parameter 

H = hu/E, (16) 

whose value characterises the relative sizes of quantum and classical effects in 
what follows. In the 'Heisenberg picture' of classical dynamics, the equation 
of motion for a(t) is 



(17) 



and, since H is a constant of the motion, 



a(t) = a(0)exp(-iutH'(H )) . (18) 

If there is a degree of uncertainty in the initial conditions of the system, 
described by a Liouville density p(a,a), then the dynamics of the system 
may be more conveniently described in the 'Schrodinger picture' of classical 
dynamics, where the state at time t is described by the new Liouville function 
p(a,a, t) = p(a(—t),a(—t)) [T7j. When K > 1, so that if is a nonlinear 
function of Ho, the dynamics leads to rotations in the complex a, a plane 
that vary in angular frequency as a function of the radial distance from 
the origin. The effect of this radial variation in the angular frequency on 
an initial density profile localised near some point in the phase plane, is to 
create 'whorls' about the origin ^H] as time progresses, as shown in Fig. ^ 
for the case H = Hi/ E. 

We are interested in the statistical properties of the evolving classical 
state, such as the mean values of position and momentum and the variance in 
these quantities. In general these statistics must be determined numerically, 
but the calculations may be simplified for certain classes of Liouville densities. 
We consider Gaussian densities 

P7,(9o,po)(^p) = 7^ ex P(-(! m ^ 2 (<?-<?o) 2 + \(p ~ Pof / m) h 2 ) , 

which can be rewritten as 

P K ,a (oi) = —exp(-K,\a - a \ 2 ) , (19) 

where k = fku/'y 2 , «o = (y/m^qo + Wo I y/TcuS) / \f2h, and the effect of the 
change of measure has been incorporated into the normalisation co-efficient. 

When the initial density is of this form, the mth moment of a at time t 
is given by 

(a m )(t) = - [ a m e- iujmtH ' { ^ a) e- Kla - aol2 d 2 a. (20) 

This integral can be simplified by setting a = rexp(z0)/y / 2 and «o — 
r o exp(20 o )/A/2 and then evaluating the integral over so that 

(a m ){t) = e ^o—g/2 / r m+i e -Kr V* Im ^rro)e- imujtH ' {hwr2l2) dr , 
2 m l J 

(21) 

where I m denotes the modified Bessel function [XHJ of the first kind of order 
m. The classical moments of primary interest are the mean and standard 




Figure 1: (Color online) Density plots showing the classical evolution of an 
initial Gaussian density centered at a® = q$ = 0.5, with k — 2, as generated 
by the Hamiltonian H = Hq/E. The parameters m, u, E have been set equal 
to 1, and the times of the plots are, from left to right and top to bottom, 
t = vr/4, t = vr/2, t = 3vr/4 and t = vr. 



deviation in position and momentum. These can be constructed from the 
moments in a and a as 



(q)(t) = y/ (2h/muj) Re{(a) (t)} , (p)(t) = V2hmu Im{(a)(t)}, 
Aq(t) = vto«|a| 2 > - Re{(a)(t) 2 })/mu - (q)(t) 2 , 



Note that the moment (|a| 2 )(t) = |«o| 2 + 1/k is a constant of the motion. 

If K > 1, then it follows from (J2U) that (a m )(t) -> as t -> oo for 
every m > 0, by the Riemann-Lebesgue Lemma ^Hl- in particular, the mean 
position and momentum tend towards zero as time increases. This is in sharp 
contrast to what happens with the corresponding quantum evolution, where 
periodic behaviour occurs 

On phase space, the quantum dynamics is determined by the Moyal 
bracket expansion (fTTj) , from which comparisons can be made with the classi- 
cal dynamics in the phase space setting. As already indicated, the transition 
from classical to quantum dynamics can be studied by successively adding 
on to the classical Poisson bracket evolution, higher order terms in that 
expansion 0, [3] until the full quantum dynamics is obtained. Note that 
for polynomial Hamiltonians such as (|14j). the series terminates. This 
process then defines a terminating sequence of semiclassical approximations, 
starting with the classical evolution, and ending with the quantum one. 

On the other hand, in order to define semiquantum approximations to 
classical dynamics, we represent both the quantum and classical dynamics 
on Hilbert space in terms of the Hamiltonian operator H which, given (|14j). 
takes the form 



The dimensionless co-efficients are determined by, but are not identical to, 
the bk in (|14jh except that ck = bj<. For example if H = Hq/E 2 , then H = 
HQ/E 2 +5(hu) 2 Ho/4:E 2 as in (|4()|) below. The relation between the coefficients 
bk and q can be determined using recurrence relations, but explicit formulas 
are very complicated. What is important to note here is that, with H a 
'classical' Hamiltonian and the bk assumed independent of h, the Cfc for k > 
typically have the form = b^ + o(h) as h — * 0. 

Since H is a function of the oscillator Hamiltonian operator Hq, it can 
be diagonalised on the well-known number eigenstates. We introduce the 



Ap(t) = y/hmu((\a\ 2 ) - Re{(a)(t) 2 }) - (p)(t) 2 . 





(23) 



k=0 



creation and annihilation operators a = W and a) = W and the 



on the states \n), n = 0, 1, 2, . . . with TV |n) = n\n), and has eigenvalues 
given by 



Our object is to construct semiquantum approximations by using the 
expansion (JZj) and, beginning with the quantum dynamics, to successively 
add on higher order correction terms to this until the full classical dynamics is 
obtained. Because the series (J7J) also terminates for polynomial Hamiltonians, 
we get in this way a terminating sequence of semiquantum approximations, 
beginning with the quantum evolution, and ending with the classical one. 

4 Example 1: Hamiltonian of degree 4 

The first example that we consider has 



H = H$/E, H = Hl/E + (fku) 2 /4 = fihw{N 2 + N+ 1/2) . (25) 



Since H is quadratic in H 0} and hence quartic in q and p, the dynamics is 
nonlinear. It is then to be expected that the classical and quantum evolu- 
tions produce comparable values and behaviours of the expectation values 
of observables for only a limited period of time, sometimes referred to as 
the break time [2D]- This depends on h and the properties of the initial 
state. The numerical value attributed to the break time in a particular case 
is dependent on which criterion is used to determine differences between the 
quantum and classical dynamics, so it should be interpreted only as a guide 
to the time-scale over the which the evolutions produce similar values for 
observable quantities. 

Differences between the classical and quantum evolutions are already ap- 
parent if we compare the classical picture in Fig^ with the quantum evo- 
lution of the Wigner function on phase space in FigJ2l which starts with 
the same initial Gaussian density, and covers the same length of time. In 
the first place, the Wigner function immediately develops negative values on 
some regions, shown in white in FigO Secondly, the quantum evolution is 
periodic, unlike the classical evolution, with period 27r//xo> in this case. This 
can be seen more clearly from the formula (|27jl below for the matrix elements 
of the density operator, all of which have this period. Many studies in recent 
years have explored in detail such characteristic differences between quantum 




Then H is diagonal 




(24) 



Figure 2: (Color online) Quantum evolution of an initial Gaussian Wigner 
function, with the same parameter values used in Figure ^ and shown at 
the same times. Regions on which the pseudo-density becomes negative are 
shown in white. 

and classical dynamics, especially in the context of 'chaos' ; see for example 

I2DE2]. 

In this first example, the full classical evolution on Hilbert space is given 
by adding just the first correction to the quantum evolution in J7J), and hence 
in this case there are only the zeroth order semiquantum dynamics and the 
full classical dynamics to study. Similarly, the full quantum evolution on 
phase space is given with just the first correction to the classical evolution 
in (jllj) . so there are only zeroth order semiclassical dynamics and the full 
quantum dynamics to study. In short, there is no dynamics 'in between' the 
classical and quantum dynamics in this case. Nevertheless, the simplicity of 
this example makes it relatively easy to derive analytic expressions for the 
classical dynamics on Hilbert space that are relevant to more complicated 
cases. 

The quantum evolution in Hilbert space 7i in this example is determined 



by 

ft = ^[^ + l),p}^-ipu[N(N + l),p]. (26) 

If we represent p in the number basis, its matrix elements p nm = (n\p\m) 
then evolve as 

Pnmit) = 6^-^(0) , (27) 

where E n , E m are eigenvalues of H with, from (J23jl. i? n = pftw(n 2 + n+ 1/2). 

Changing variables from (g,p) to (a, a), we find the classical evolution 
(jlj), on 7i in this example as 

— = -^[iV(iV + l),G] 

where = W^H/da 1 ), H a7 , = W" 1 (d 2 H / dada) , etc. In this expres- 
sion the correction term appears to be of lower order in h than the zeroth 
order term because we have introduced the operator N which is 0(1/ h), and 
derivatives with respect to a and a which are 0(y/~H). After some algebraic 
manipulation, (|2*8|) simplifies to 



1 

— = -i(puj/2) (\N(N + 1), G) + [N, aGa) + a f Ga]) . (29) 

Note that although as remarked above this must be fully equivalent to the 
classical evolution of the Liouville density, there is still dependence on h in 
the RHS, because of the dependence implicit in the definition of G. The 
effect of the evolution (J29|) can be made clearer if we write G as the 'operator 
vector' \G)) G TL (&TC* jH] and replace the right and left action of operators 
on G by superoperators [2311211113 acting on \G)), labelled by the subscripts 
r and I respectively, i.e. A{\G)) = AG and A r \G)) = GA^. We can then 
rewrite (j29| in the operator vector form 

-i(puj/2)(N l - N r ) (Ni + N r + 1 + aia r + ajaty \G)) . (30) 



d\G)) 



dt 

If we now define the new superoperators 

AL = N t - N r , X x = (N t + N r + l)/2 



X 2 = (d[a r + a\a\)/2 , X 3 = (aia r — a\a\)/2 , (31) 



then the close to form an sl(2) algebra of superoperators that commute 
with N_ and satisfy on TC ® TC* the commutation relations 

[X 2 , XJ = X 3 , [1 3 , X 2 ] = X, , [X 3 , X,] = X 2 . (32) 

In terms of these new superoperators, the classical evolution of \G)) is given 
by 

d\G))/dt = -iiiuN_(X 1 + X 2 )\G)) . (33) 

Since the sl(2) superoperators commute with iV_ , one method of simplify- 
ing the evolution is to decompose the tensor product space TC®TC* into repre- 
sentations of si (2) labelled by the eigenvalues of N_ . These eigenvalues may 
take any integer value ^GZ and the corresponding si (2) lowest- weight oper- 
ator eigenvector of X\ is given by \u, 0)) for v > and |0, —v)) for v < 0. Here 
we have introduced the basis of operator vectors {\n, m)), n, m = 0, 1, 2, . . .} 
in TC <g> that corresponds to the basis of operators {|n)(m|, n,m — 
0, 1, 2, . . .} in the space of operators acting on TC, with A^|n, m)) = n\n, m)) 
and N T \n,m)) = m\n, m)). 

Now we can write 

(34) 



v > 

(35) 

1/ < 
1/ > 0, 

G^ = ((n,n- v\G)) = (n\G\n - v) for v < . (36) 
The classical evolution (|3*3*j) of \G)) generated by H now leads to 

\G{t))) = £|G"(f)» = ^exp(-z^tpM)|GM(0))), (37) 

fGZ i/GZ 

where the only nonzero matrix elements of the superoperator p( v > are given 
by 

(P») nn = n+(|i/| + l)/2, 
(P M )„„+i = (P M )„+i„ = y/(n + l)(n+\u\ + l) 12. (38) 



|G» = JJ |GM» , 
where 

|GM» = 

I Er=0 G nV>™-^)> » 

with 

= (( n + u,n\G)) = (n + u\G\n) for 



As a counterpoint to the quantum evolution ()27|) of the matrix elements of 
p, we find that the classical evolution transforms the matrix elements of G 
along each diagonal independently: 

oo 

G nm {t) = V (exp[-i(n - m)//u;tP (n - m) ]) G„„ m+r . r (0) , (39) 

r=0 

for n > m. By the z4h diagonal, we mean the set of matrix elements for which 
the row label minus the column label is equal to v. If G is initially Hermitian, 
then it remains Hermitian under this evolution and the matrix elements of 
G(t) above the main diagonal can be obtained by complex conjugation from 
those below. 

We can deduce at once several properties of the classically evolved matrix 
G(t) by recalling some properties of the classical evolution on the phase plane. 
Firstly, the total amount of probability on the plane remains constant, which 
implies that the trace of G(t) should also remain constant. That this is 
reproduced in the Hilbert space analysis can be seen from ()39p by setting 
n = m and observing that the diagonal elements of G are constants of the 
motion. Secondly, the integral over phase space of the square of the Liouville 
density is also a constant of the motion, and this is equivalent to the Hilbert 
space condition that Tr(G(t) 2 ) should remain constant. This property of 
G(t) can be verified from (|39|) by noting that the evolution of each diagonal 
of G(t) is unitary and hence that the sum of the squares of the elements of 
the z/th diagonal is a constant of the motion. It then follows that the sum 
over all diagonals, equal to Tr(G(t) 2 ), is also a constant. 

Note however that the classical evolution generated by H is not unitary 
in the sense one uses when describing quantum dynamics, that is to say, in 
the Hilbert space 7i. Instead, it defines an evolution of G that corresponds to 
a unitary superoperator evolution of the operator vector \p)) in 7i®7i* . This 
has the consequence that some features of a quantum evolution remain, such 
as the trace relations described above, but that other characteristic features 
of a quantum evolution are lost. For example, there is no requirement for the 
eigenvalues of the Groenewold operator to remain constant under the motion. 
These eigenvalues can and do change and in general, negative eigenvalues 
will typically develop over time, even if all eigenvalues of G are intially non- 
negative. 

When examining the relationship between quantum and classical dynam- 
ics, it is primarily of interest to consider the limit h — > in some appropriate 
way. There is some flexibility in how to do this. In ^H], the limit was taken in 
such a way that the initial phase space density approaches a delta-function in 
q and p , so that both the quantum and classical dynamics approach that of 



a classical trajectory. In this paper, we shall instead investigate numerically 
the effect of reducing the value of h in such a way that the initial Gaussian 
phase space density is kept unchanged. 

This involves a scaling of three key parameters: re, \x and olq. The first 
two of these depend linearly on h, while «o is linear in (1/y/h). In order 
to keep the initial density constant as ft — > 0, we must have re, fi — > and 
a — > oo, while keeping re/// and re|a;o| 2 constant. 

In FigOl we graph the first and second moments of a for two sets of the 
parameters re, \x and ao, with t in the range [0, 7r]. For the top two graphs 
in FigOl the values re = 2, (/, = 1/2 and ao = 1/2 are used, while in the 
bottom two graphs, the values re = = 1/4 and a = l/y/(2) are used, 
corresponding to H — > H/2. The first moments are graphed on the left in FigC2 
and it can be seen that the quantum evolution lies close to the classical for 
much longer with the smaller (effective) value of h. A more exaggerated sign 
of this convergence can be seen in the graphs of the second moments on the 
right. In the upper graph, one can see that the classical moment approaches 
the origin, whereas the quantum moment exhibits a symmetry in q. This 
symmetry is a sign of the recurrence of the initial state at t = 2n. In the 
lower graph however, the recurrence is not due to occur until t = 4n and the 
quantum evolution resembles the classical much more closely for t G [0, 2tt]. 

5 Example 2: Hamiltonian of degree 6 

In order to compare and contrast semiquantum and semiclassical dynam- 
ics, we need to select a Hamiltonian of least degree 3 in Hq, so that there 
are nontrivial semiclassical and semiquantum approximations lying between 
classical and quantum dynamics. The simplest such case has 

H = Hl/E 2 , H = Hl/E 2 + 5{hu) 2 H /4E 2 

= j2 2 huj(N 3 + 3N 2 /2 + 2N + 3/4). (40) 

The eigenstates of H are again the number states and it is natural to 
represent the dynamical evolution in Hilbert space in this basis. The quan- 
tum evolution of the matrix elements of an initial density operator p(0) on 
7i is again given by (|27|h where in this case the energy eigenvalues take the 
form E n = fi 2 hw(n 3 + 3n 2 /2 + 2n + 3/4). From this we see that the quantum 
evolution is again periodic, but now with period 47r/ fi 2 u. 

The equations governing the classical evolution of the Groenewold oper- 
ator on 7i are obtained as before by inserting the Hamiltonian into equation 




Figure 3: (Color online) Classical and quantum evolution of the first and 
second moments of a = (q + ip) / \/2 as generated by H = Hq / E. Points 
on the classical curves are marked by + and points on the quantum curves 
by x. Again m = uj = E = 1, and here go = V%, Po = 0. The first 
moments are graphed on the left and the second moments on the right. In 
the two upper graphs, k — 2, /x = 1/2, a = 1/2 and in the lower two graphs, 
« = l, Ai = l/4,a = l/ v ^2). 



(|7|). The result is that we obtain an expression for the full classical evolution 
that starts with the quantum evolution and adds two correction terms: 



dC r 

— = -ip 2 hw N 3 + 3N 2 /2 + 2N, G 
dt L 



Note that terms involving fourth derivatives with respect to a and a do not 
appear in the second classical correction, because <9 4 if/<9a 4 and d 4L H/da i 
vanish in this example. 

The correction terms can be simplified by straightforward but lengthy 
algebraic manipulation. Of particular interest is the semiquantum evolution 
equation that includes the quantum evolution and the first classical correc- 
tion. This equation simplifies to 

dC / 

-i(/i 2 Huj/A) [5[N 2 , a)Ga] + 5[N(N + 2), aGa)} 



dt 

-[N(N + 1)(N + 2),p] - [N,a 2 pa j2 + a t2 pa 2 + 3NpN]^j . (42) 

If we make use of the si (2) operators (|3*T|) introduced earlier, then the 
corresponding evolution equation for the operator vector \G)) is given in this 
case by 

d\G))/dt = -i^ujN. (3(XxX 2 + X 2 X x )j2 - (X t - X 2 ) 2 ) \G)} . (43) 

The solution of this equation has similar properties to the solution of (|HBJ). 
In particular, it can again be decomposed as in (JBH), leading in this case to 
an expression of the form 

\G(t)}} = £|G"(t))> = J2 G M-i^tQ iu) )\G v (0))) (44) 

in place of (|37|) . However in this case, truncation methods for determining 
the eigenvalues of Q v appear to fail, and hence for computational purposes, it 
is useful first to introduce the unitary si (2) operator U = exp(log(7/3)Xs/4) 
which enables us to simplify the expression in ()43j) . as 



U (3(XiX 2 + X 2 Xi)/2 - {X, - X 2 ) 2 ) fr = V2i [x x X 2 + X 2 xA . (45) 



This is seen after noting from (J3T|) and that 

3(XxX 2 + X 2 X x )/2- (Xi- X 2 ) 2 = 3X1/4 -7X1/4, 

X 1 X 2 + X 2 X 1 = (X 2 + -X 2 )/4, (46) 

where 

X± = X 2 ±X 1 , X 3 X ± = X ± (X 3 ±1). (47) 

The transformation (|43j) maps to a tri-diagonal matrix QW, whose only 
non-zero matrix elements are given by 

Q { Sn = Qnn + i = V21(n + 1 + \u\/2W(n + l)(n+\u\ + l)/2, (48) 

and for which the eigenvalues can easily be determined numerically. The 
expression for the individual matrix elements of the Groenewold operator in 
this semiquantum approximation is now given by 

oo 

G nm (t) = U t~ m) exp (-i(n - m)^utQ^') U%- m *G n - m+t M , 

r,s,t=0 

(49) 

where Unm = (n + u, n\U\m + u, m) for v > 0. 

In order to determine the full classical evolution on Hilbert space, we 
need to evaluate both correction terms in (j41j) . This leads in place of (jHSJ) 
to a comparatively simple expression for the evolution of the operator vector 

d \ G )) o,..2 /AA y t * , 



dt 



^ l //AN.{X 1 + X 2 ) 2 \G)) . (50) 



Note that this evolution again depends on the operator sum X\ + X 2 . 
This property turns out to be shared by the classical evolution for each 
Hamiltonian of the form (JHJ). The formula corresponding to (|5Uj) in the 
general case is 

m. = -i^N. ± (*±^(*. + X^G)) . (51) 

k=0 



However, expressions for the successive semiquantum approximations to the 
quantum evolution, such as (jEJ), are not so easy to obtain for a general 
Hamiltonian of the form (fT4*|) . 



Returning to the example at hand, with H as in ([40)1. the classical evo- 
lution of the matrix elements of G is given by 

G nm (t) = (ex P H(n - m)^ut (p>-™)) ] J G n _ m+r , r (0) , (52) 

r= Q \ / mr 

with P^) as in (|38j). Eqn. (|52|) is to be compared with (|39|) in the previous 
section. Working from (j4*9|l and (152(1 . we can evaluate and compare various 
quantities of interest in the semiquantum and classical evolutions, such as 
the expectaton values of a, and the eigenvalues of the Groenewold operator. 

There is also a semiclassical approximation to be studied in this case. The 
classical dynamics involves two correction terms to the quantum evolution, 
and the quantum dynamics involves two correction terms to the classical 
evolution. Just as the semiquantum dynamics includes the first correction 
to the quantum dynamics but not the second, so the semiclassical dynamics 
includes the first correction to the classical dynamics but not the second. It 
is of particular interest to determine differences between the semiquantum 
and the semiclassical dynamics. 

In Fig. (|^), we show phase space plots of the evolution of the expectation 
value of a in the semiclassical and semiquantum approximations, along with 
the classical and quantum evolutions. The parameter values were chosen to 
be k = 2, fj, = 1/2 and Oq = 1/2, and the time evolution is over the interval 
[0, it]. It is immediately clear that the semiclassical and semiquantum approx- 
imations are quite distinct. For small values of t, they closely approximate 
the quantum and classical moments respectively. However, the long-term be- 
haviour is very different, as the semiquantum moment begins to oscillate in 
a fashion that is qualitatively similar to the quantum moment, whereas the 
semiclassical moment appears to approach smoothly a fixed constant, just as 
is the case for the classical moment. 

How does the behaviour change as H decreases? To examine this, we 
graph in Fig.© the first moment in a over t G [0, tt], but this time with 
the parameter choices k = l,/i = 1/4, a — l/v2~, which is equivalent to 
dividing h by two. The first thing one observes is that the quantum evolution 
is now much closer to the classical, although it begins to oscillate near t = tt. 
It is also apparent that the semiquantum and semiclassical curves are good 
approximations to the classical and quantum moments respectively over most 
of the time-interval. One can also add that the semiclassical moment appears 
to be a good approximation to the quantum moment for a longer time than 
the semiquantum moment is a good approximation to the classical moment. 

All the forms of dynamics can be studied in phase space or in Hilbert 
space. However, it must be emphasized that here as in general, the zeroth 
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Figure 4: (Color online) Comparison of first moments of a = (q + ip) / \[2 for 
classical, semiquantum, quantum and semiclassical evolutions generated by 
the Hamiltonian H = Hq/E 2 . Points on the classical, quantum, semiclassical 
and semiquantum curves are labelled by +, x, o and □, respectively. The 
evolution is over the time-interval [0, it] and again m = oj = E = 1, with 
k = 2, fi = 1/2, ceo = 0.5. 
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Figure 5: (Color online) Identical to the preceding Figure except that h has 
been reduced by a factor of two, by setting k — 1, jx — 1/4, olq — l/y/2. 



order semiclassical approximation to quantum dynamics, namely classical 
dynamics itself, when applied in Hilbert space with an initial Groenewold 
operator G(0) that is the image W _1 (27r/ip) of an initial Liouvillean density, 
will in general produce a G(t) that is not positive-definite for times t > 0, even 
if G(0) is positive-definite, which it need not be. Closely related observations 
have been made in earlier studies of the evolution of density operators under 
classical dynamics 1231 121 123 . 

Conversely, the zeroth order semiquantum dynamics, namely quantum 
dynamics itself, when applied with an initial Wigner function W(q, p, 0) 
that is the image W(p/27rft) of an initial quantum density operator, will in 
general produce a W(q,p,t) that is not positive-definite for times t > 0, 
even if W(q, p, 0) is positive-definite, which it need not be. 

The preceding figures use measures of comparison that are natural in the 
phase space formulation of all the dynamics. Alternatively, one can use mea- 
sures that are typically associated with formulations on Hilbert space. In 
particular, one can examine the evolution of the eigenvalues of the density 
operator (or the Groenewold operator) in the classical, quantum, semiclas- 
sical and semiquantum approximations. In Fig.© we show the evolution of 
the largest pair and smallest pair of eigenvalues for the classical, semiquan- 
tum and semiclassical cases. In the quantum case, the eigenvalues are and 
1 at all times as the system is in a pure (coherent) state initially, and stays 
in a pure state at all subsequent times. The parameter choices and time 
interval used in Fig.© correspond to those used in Fig.©. We also single 
out the points at which t = 1, 2, 3 on each graph for the purpose of compar- 
ison with the moment curve of Fig.(j3J) which is reproduced in Fig.©. From 
these graphs, we observe that the semiquantum approximation to the classi- 
cal moment is reasonable up till about t = 1, by which time the eigenvalues 
of the classical and semiquantum Groenewold operators differ to a significant 
degree. The semiquantum eigenvalues appear to display periodic behaviour. 

The semiclassical eigenvalues remain very close to the classical eigenvalues 
over the whole time range. Since the behaviour of the individual eigenvalues 
does not seem to reflect the differences shown in the graphs of the first 
moment, this leads one to wonder if those differences are reflected in some 
global property of the spectrum. We know that in each case, the evolved 
Groenewold operators have trace and square trace equal to 1, but we can look 
at the contribution to these quantities by the negative eigenvalues. Numerical 
experiments indicate that it is difficult to calculate the sum of the negative 
eigenvalues by using truncation techniques, so we instead concentrate on the 
sum of the squares of the negative eigenvalues. In Fig.© we graph this 
'squared negativity' for the parameter values used in Fig.fjlJ) 

These results, however, merely confirm the conclusions drawn from the 
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Figure 6: (Color online) Comparison of largest two and least two eigenvalues 
for, from left to right and top to bottom, classical, semiquantum and semi- 
classical evolutions generated by H — H^/E 2 , for the time interval [0, 7r], 
and with m — oj — E — 1 and n = 2,fi = 1/2, a = 0.5. The evolution of 
the first moment is reproduced from Fig|3] in the graph at bottom right for 
comparison. Each of the other graphs also features the quantum spectrum 
{0, 1} and in all graphs the values at the time-points t = 1, 2, 3 are marked 
(classical), o (semiclassical), A (quantum) and □ (semiquantum). 
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Figure 7: (Color online) Comparison of 'square-negativity' for classical, semi- 
quantum and semi-classical evolutions generated by H = Hq/E 2 for the 
time-interval [0, w], and with m = u = E = 1 and k = 2, fi = 1/2, a = 0.5. 
The solid line, the solid line with x points and the solid line with + points 
represent the classical, semiquantum and semiclassical values respectively. 



eigenvalue curves and indicate that the spectrum of the semiclassical Groe- 
newold operator resembles that of the classical Groenewold operator more 
closely than the spectrum of the semiquantum Groenewold operator does, for 
all but very small times. We conclude that the moments and the eigenvalue 
spectra provide quite different information about these different approximate 
dynamics. 

6 Concluding remarks 

Our study indicates that semiquantum approximations to classical dynamics 
can provide new information about the interface between quantum and clas- 
sical mechanics. Semiquantum approximations obtained from the Hilbert 
space formulation of classical mechanics, and semiclassical ones obtained 
from the phase space formulation of quantum mechanics, show significantly 
different behaviours for expectation values of observables and eigenvalues of 
the (pseudo) density operator. For the simple systems considered here, it 
appears that a first semiquantum approximation to quantum dynamics is, 
for simple indicators like moments of coordinates and momenta, closer to 
classical dynamics and further from quantum dynamics than is a first semi- 
classical approximation to classical dynamics. On the other hand, as regards 
the spectrum of the density operator and Groenewold operator, the first semi- 
classical approximation behaves more like classical dynamics than the first 
semiquantum approximation, which behaves more like quantum dynamics in 
this respect. It is not at all clear why this should be so, and it is desirable 
in future work to look at cases where higher approximations come into play 
to see what happens then, as well as to examine the underlying theory more 
closely. 

We have been able to consider only very simple one-dimensional systems 
here, and there is obviously a need also to explore systems with more degrees 
of freedom, especially nonintegrable ones. The possible role of semiquantum 
mechanics in throwing new light on the interface between classical chaos and 
corresponding quantum dynamics is especially interesting. 

More generally, there are deep theoretical questions that arise about the 
mathematical relationship between semiquantum and semiclassical approxi- 
mations, associated with the structure of the series expansions (J2J) and (jlljl . 
Another concerns the nature of the relationship between action principles and 
semiquantum approximations to classical dynamics, whether formulated in 
phase space or Hilbert space. We hope to return to some of those questions. 
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